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O ' Abstract 

We present an algorithm to solve First Order Ordinary Differential Equations 

^.^ (FOODEs) extending the Prelle-Singer (PS) Method. The usual PS-approach miss 

r~^ ■ many FOODEs presenting Liouvillian functions in the solution (LFOODEs). We 

^^ ! point out why and propose an algorithm to solve a large class of these previously 

^p I unsolved LFOODEs. Although our algorithm does not cover all the LFOODEs, it is 

f~^ ■ an elegant extension mantaining the semi-decision nature of the usual PS-Method. 
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1 Introduction 

The problem of solving ordinary differential equations (ODEs) has led, over the 
years, to a wide range of different methods for their solution. Along with the 
many techniques for calculating tricky integrals, these often occupy a large part of 
the mathematics syllabuses of university courses in applied mathematics round the 
world. 

The overwhelming majority of these methods are based on classification of the 
DE into types for which a method of solution is known, which has resulted in a gamut 
of methods that deal with specific classes of DEs. This scene changed somewhat 
at the end of the 19th century when Sophus Lie developed a general method to 
solve (or at least reduce the order of) ordinary differential equations (ODEs) given 
their symmetry transformations [Q, ||, |3[. Lie's method is very powerful and highly 
general, but first requires that we find the symmetries of the differential equation, 
which may not be easy to do. Search methods have been developed |Q, ^ to extract 
the symmetries of a given ODE, however these methods are heuristic and cannot 
guarantee that, if symmetries exist, they will be found. A big step forward in 
constructing an algorithm for solving first order ODEs (FOODEs) analytically was 
taken in a seminal paper by Prelle and Singer (PS) |^ on autonomous systems of 
ODEs. Prelle and Singer's problem is equivalent to asking when a FOODE of the 
form y' = M{x,y)/N(x,y), with M and A^ polynomials in their arguments, has 
an elementary solution (a solution which can be written in terms of a combination 
of polynomials, logarithms, exponentials and radicals). Prelle and Singer were not 
exactly able to construct an algorithm for solving their problem, since they were 
not able to define a degree bound for the polynomials which might enter into the 
solution. Though this is important from a theoretical point of view, any pratical 
use of the PS method will have a degree bound imposed by the time necessary to 
perform the actual calculation needed to handle the ever- more complex equations. 
With this in mind it is possible to say that Prelle and Singer's original method is 
almost an algorithm, awaiting a theoretical degree bound to turn it algorithmic. 
In 0], we have presented an extension of the PS-ideas to deal with second order 
ordinary differential equations (SOODEs). The purpose of this paper is to present 
another extension to the Prelle-Singer method allowing for the solution of a class 
of LFOODEs. The paper is organized as follows: in section ||, we present a short 
theoretical introduction to the PS approach; in the following section, we show why 
the Prelle-Singer aproach misses some LFOODEs; in section ^ we introduce the 
main ideas of our extension with examples of its application; we finally present our 
conclusions. 

2 The Prelle-Singer Procedure 

Despite its usefulness in solving FOODEs, the Prelle-Singer procedure is not very 
well known outside mathematical circles, and so we present a brief overview of the 
main ideas of the procedure. Consider the class of FOODEs which can be written 

as 

, ^dy_ ^ M{x,y) 

^ dx N{x,y) ^ ' 

where M(x,y) and N{x,y) are polynomials with coefficients in the complex field 
C. If i? is a integrating factor of {^, it is possible to write dx{RN) + dy{RM) = 0, 



leading to N d^R + R dyN + M dyR + R dyM = 0. Thus, we finally obtain: 

^ = - {d.N + dyM) , (2) 

where D = Ndx + Mdy. In |^], Prelle and Singer showed that, if the solution of ([l|) 
is written in terms of elementary functions, then R must be of the form -R = Hi /r"' 
where fi are irreducible polynomials and rii are non-zero rational numbers. Using 
this result in (0), we have 

D[R] _ Diujn _^^f^'^^D[fm,^J;^ 



R n./r Ukfk' 

From (|2|), plus the fact that M and N are polynomials, we conclude that D\R\/R 
is a polynomial. One can then prove from (|3[) that fi\D[fi] ^. We now have a 
criterion for choosing the possible fi (build all the possible divisors of D[fi\) and, 
by using ^ and (y), we have 

Y."^^ = - {d.N + dyM) . (4) 

If we manage to solve @ and thereby find ni, we know the integrating factor for 
the FOODE and the problem is reduced to a quadrature. 



3 The PS-Method and LFOODEs 

In the previous section, the main ideas and concepts used in the Prelle-Singer 
method were summarized. In this section, we will first explain why the PS-method 
fails to solve some LFOODEs. Then, using the results on Q , we will show that it 
is possible to algoritimicaly solve a large class of these. 

3.1 Introduction 

As we have already mentioned, the usual PS-method garantees to find a solution 
for a FOODE if it is expressible in terms of elementary functions. However, the 
method can also solve some FOODEs with non-elementary solutions. Why and 
why not for all the FOODEs with non-elementary solutions? Consider the following 
two examples: 

dy 3x^y^ + x^ + l 

dx 4(x + 1) (x2 -x + l)y 
and 

y +yx + x-l, (6) 



dx 

equation 1.18 of the standard testing ground for ODE solvers by Kamke |^]. These 
equations present (respectively) general solutions given by: 



Va; + l\/x2-x+l (l/2 j- 



y^ - V^TTVx2 -x+l ( 1/2 / }^— dx - 2_ci ] = 0, (7) 

\/x + \\/x^ — x + l / 



and 

_ 2_C1 + ./^^e-^^/2erf{l/2^J^^/2x - ^/^V2) - 2e^/^^(^-^) 
^~ 2 _C1 + V^\/vFe-2v^er/(l/2 ^/^^/2x - ^/^\/2) ' ^^ 

We note that both solutions are not expressible in terms of elementary functions. 
But, for F00DE|5|, the standard PS-method can find the solution (eq. 0). The same 
is not true for FOODE |6|, what is happening? This can be best understood if we 
have a look on the integrating factors for those FOODEs, which are respectively: 

i?= (x3 + l)"^^^ (9) 

and 

z2/2-2a; 

(y + 1) 

Since the standard PS procedure constructs integrating factor candidates from poly- 
nomials in the variables {x,y), one can see that, since the integrating factor on eq. 
IC presents the exponential e^ I'i-'ix ^ j^. ^-^ never be found by the PS-method. So, 



can we understand something about the general structure of the FOODE that will 
allow us, eventually, to solve a class of these latter equations? 

3.2 Important results concerning the general form of 
the integrating factor for LFOODEs 

In this section, we are going to present a result concerning the general form for the 
integrating factor for LFOODEs. Based on that, we are going to develop, on the 
next section, a method of solving a class of those LFOODEs. Using the following 
theorem, which has been presented on |^], Theorem 1: Consider a LFOODE of 

the form dy/dx = M{x,y)/N{x,y), where M(x,y) and N{x,y) are polynomials 
in {x,y). Its integrating factor must be of the form e^J ^ ^' , where U and 

V are rational functions on {x,y) with Uy = Vx- we can extract following result: 

Theorem 2: Consider a LFOODE of the form dy/dx = M{x,y)/N{x,y), where 

M{x,y) and N{x,y) are polynomials in {x,y). Its integrating factor must be of the 
form: 

n 

R = e'''^-^y^Y[p,ix,yr^. (11) 

i=l 

where tq is a rational function on {x,y), the pi's are irreducible polynomials on 
{x,y) and the Cj's are constants. 



Proof: From well known results concerning formal integration |1C], we have 
that if f{x) is a rational function of x, we can conclude that: 



K 
f{x)dx = ro{x)+Y,aHr,{x)) (12) 



where rj {i = 0, ..., K) are rational functions of x, and C, (i = 1, ..., K) are constants. 
Let U{x,y) and V{x,y) be two functions such that Uy = Vx- So, we can choose 



a function uj(x,y) such that dto = Uxdx + Vydy, i.e., we can write uj^ = U and 
LUy = V. Integrating this system (to obtain uj) we have 



CO 



fu dx+ fvdy- [ dy ju dx dy. (13) 



The integrals appearing above are what we call partial integrals (otherwise, equation 
(|l3 ) would not be true). Basicaly, a partial integral is performed by considering all 
but the integration variable as constants. Therefore, if we also use that U{x, y) and 
V{x, y) are rational functions, one can use the result ( [l2|) thus obtaining: 



jUdx = UQ{x,y)+Y,C^Hu,{x,y)) (14) 



and 

K 



Jvdy = vo{x,y) + Y,C^ln{vi{x,y)) (15) 



where «i (i = 0,...,K^), vi {i = 0, ...,-fC^) are rational functions of x,y, and C^^ {i 



1, ...,K'^),Cl' {i = 1,...,K'") are constants. We also have, using (|14D : 



dy ju dx = dyuo + ^ C," ^. (16) 



It is straightforward that the righthand side of the above equation is also a rational 
function on (x,y). So, J dy JU dx dy can be written on a similar form as presented 
on eqs. ( |l4| , |l5| ). Thus, using all these results, we can conclude that w can be writen 
as: 

where Wi {i = 0, ...,K^) are rational functions of x,y, and C^" {i = 1, ...,K^) are 
constants. So, combining theorem 1 and eq. ([l7|), it is straightforward to see that: 



R = e(J(Udx+Vdy)) ^ ^ij(dujix,y))) ^ ^uj{x,y) ^ ^wo{x,y)+J2fZ ^T HMx,y)) (^g) 

resulting in: 

R = e'"<^^''^y'^Y[wi{x,yf^ (19) 

i=l 

We can write each Wi{x,y) as: 

rii 

Mx,y) = l[pj(.x,yp (20) 

where the Pj{x,y) are 'monies' in {x,y). We can write, finally: 

i? = e^«("'^)np,(x,y)^^ (21) 

Notice that, if ro{x,y) is constant, the integrating factor R would be within reach 

of the PS-method. In what follows, we are going to present a method to find R for 
a class of equations where ro(x, y) is not a constant (so the usual PS-method fails). 



4 Finding the Integrating Factor 

In this section, we introduce a method of finding the integrating factor for a class 
of LFOODEs using an algorithym somewhat similar to the PS-method. First we 
analyse the generality of our method. We then detail it and finally present some 
examples of its applicability. 

4.1 Introduction 

On the previous section, we have deduced the general form for the integrating factor 
for LFOODEs that can be written as in eq. (|^). By using R as given on ( [2l| ) onto 
eq. @, one finds: 

CiDlp^] fdN dM\ , , 

In the PS-method, the analogous to the above equation is eq. @: 



n,D\f,\ _ fdN dM\ 



^ fi \dx dyj- 



Note that the main difference among these equations is the presence of the "extra" 
term D[rQ] on eq. (|2^) . For eq. (§), using the fact that M and N are poly- 
nomials on (x,y), one can prove that fi are "eigen-polynomials" of the operator 
D, i.e., /i|i^[/j][ol- Regarding eq. (p2|), if we had that D[ro] is a polynomial, by 
the same line of reasoning, we would conclude that pi are "eigen-polynomials" of 
the operator D, i.e., pj|L)[pj]. Conjecture 1: Consider a LFOODE of the form 

dy/dx = M{x,y)/N(x,y), where M{x,y) and N(x,y) are polynomials in {x,y), 
with integrating factor given by R = e^'^^^''^'YYi^iPi{x,yY% where rg is a rational 
function on {x,y), the pi's are irreducible polynomials on {x,y) and the Ci's are 
constants. Then we have that -D[ro] is a polynomial (D = Ndx + Mdy). Assuming 

this to be true, we are going to present a method to find the integrating factor for 
a class of LFOODEs missed by the PS-method. 

4.2 The Method 

The method we are going to introduce Inow deals with three different general forms 
for ro{x,y): 

1- ro = ro(x) 

2. ro = ro(y) 

3. ro = r(x) + s{y) 

where all of the above functions are rationals. The LFOODEs that correspond to 
these cases define the class of equations solved by our method. 



4.2.1 Case ro = ro(x) 



Reminding the reader that the D operator is defined hy D = Ndx + Mdy, eq. (22) 
will then become: 

dx Pi \dx dy J ' 

First of all, since we made a hypothesis implying that the pj's are "eigenpolyno- 
mials" of the D operator, we have: D\pi] = gi Pi, where gi are polynomials called 
eigenvalues of D. We can then calculate all the pi and associated gi up to a given 
order, in the same fashion of the PS-method. We can then write eq.(p3) as: 

. drojx) _ fdN ^ dM 
leading to: 



~l7^-(i^- if )-!:-'■ <^^) 



dro{x) _ N^ + My + Y,i Cj gi .^^. 



dx N 

In order to find ro(a;) we need first to integrate eq.(p5|): 

Nx + My + Ei Ci gi 



roix) = -J ^ ' -y^--^^^^ dx + C. (26) 

Since ro is a rational function of x only, we can determine the Cj's by imposing 
that the derivative of the righthand side of eq.(2£) with respect to y be equal to 
zero and that the logarithmic terms that might arise from the integration vanish: 
If we manage to find a set of Cj's satisfying these conditions, in the order we are 
considering, we will have found the integrating factor for the LFOODE. Otherwise, 
we then increase the order of the pj's and try again, until we succeed, in the same 
manner as in the PS-method. 

4.2.2 Case ro = ro{y) 

Analogously, for this case, eq.(p2D will become: 

dro(y) ^^aDlpi] fdN dM\ 



Following a procedure similar to the one used on section 4.3.1, we can write rod/) 
as: 

f Nx+My + J2iCi gi r , ^ .OQN 

my) = -J j^ dy + C. (28) 

Since ro is a rational function of y only, we can determine the q's by imposing that 
the derivative of the righthand side of eq.([2^) with respect to x be equal to zero and 
that the logarithmic terms that might arise from the integration vanish: In the same 
way, if we manage to find a set of Cj's satisfying these conditions, in the order we are 
considering, we will have found the integrating factor for the LFOODE. Otherwise, 
we then increase the order of the pi's and try again, until we succeed, in the same 
manner as in the PS-method. 



4.2.3 Case ro = r{x) + s{y) 
For this case, eq.(p2[) becomes: 



dx dy ^—^ Pi \ dx dy J 

Dividing the equation above by N and isolating ^^^ , one obtains: 



dx 

dr{x) Mds{y) N^ + AL + j:^ a yi 



dx N dy N 

Integrating both sides of the equation, with respect to x: 



(30) 



.(.) = -^/fdx-/^- + \+^----^-dx. (31) 



We can then solve (BTI) to find: 



N^ + My + Ei Ci yi 



ds{y) J N 



dx — r{x) 



dy f^dx 



(32) 



Analogous procedure leads to: 



dr{x) _j M —"^y- ^yy> 

dx ~ f N 



f N 



(33) 



The Cj's remain to be determined. Once again, we still have to imposed that, by 
construction, r(x) and dr{x)/dx do not depend on y and, analogously, s{y) and 
ds{y)/dy do not depend on x. By imposing that, and that both r{x) and s{y) are 
rational functions, we may hope to find a set of suitable values for Cj. If that is the 
case, we would have found r{x) and s{y) and, consequently, the integrating factor. 
If this fails, we have to increase the order of the pj's and try all over again. 

4.3 Examples 

In this section, we are going to present one example of application of our method 
for each of the cases mentioned above. 

4.3.1 ro = ro(x) 

Consider the LFOODE (eq.I.129 in g]): 

{x + I) -T-y{x) + y{x) {y{x) - x) . (34) 

For this equation, up to order 1, we have that the "eigenpolynomials" (with the 
associated "eigenvalues") are: 

• Pi = y, 91 = (x-y), 

• p2 = {X+ 1), 52 = 1- 
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So, eq.(p6|) becomes: 

{{2 + ci)y - C2 + Cl) ln(x + 1) - x - ci x. (35) 

Leading to: 

ci = -2, C2 = -2, (36) 

and, consequently to: 

ro(x)=x, ^ ^ = ^2(/+ 1)2 (37) 

4.3.2 ro = roiy) 

Consider the LFOODE (eq.1.235 in |]): 

{x + I) -—y{x) + y{x) {y{x) - x) . (38) 

ax 

For this equation, up to order 1, we have that the "eigenpolynoniials" (with the 
associated "eigenvalues") are: 



(39) 
(40) 

(41) 



• pi = y, gi = -0, 

So, eg . (|26|) becomes: 


(-l-c,)ln(2/) + | 


Leading to: 


Cl = -1, 


and, consequently to: 


V 

ro{y) = T. -^ R= — 
y 



4.3.3 ro = r{x) + s{y) 

Consider the Abel LFOODE of the first kind 

d , . y{xf{y{x) + x - 1) 

d^y^^^ = X 



^yi^) = " ^^^:^^^-~-^ (42) 



For this equation, up to order 1, we have that the "eigenpolynoniials" (with the 
associated "eigenvalues") are: 

• pi = X, gi = X, 

• P2 = y, 92 = y'^ + y*x-y, 

• P3 = X + y, g3 = x-y + y^. 



+ 



So, from eq.(|33D we get: 

drjx) ^ {{-C2 - 2) x'^ + {2 C3 + 2 C2 + 6 + Cl) X - C3 - C2 - 2) ylnjy) 
dx x^ (x — 1 + ln{y)y — ln(y + x — l)y) 

((-cg - l)x^ - Cj X- 1) j/ln(j/ + x - 1) 
x^ (x — 1 + ln{y)y — ln{y + x — l)y) 

{-x^ + 2x-l)ys{y) + {c3 + 2 + ci)x^ + {-ci " cg - 2) x ^ ^ 
x^ (x — 1 + ln(y)y — ln(y + x — l)y) 
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(43) 



Leading to: 

■ ^^y + ( (~C5 - l)x^ - ci X - l) y] ln{y + x - 1) + 



dr{x) o^^ , // „ 1^ 2 



dx 



(^^x^y + ((-C2 - 2) x2 + (2 C3 + 2C2 + 6 + c, ) X - C5 - C2 - 2) y^ ln(y) + 



(x^ - x^) ^j^ + (-x^ + 2x - 1) ys(y) + 



dx 

„2 



(c5 + 2 + Ci)x^ + (-cj -C5-2)x = (44) 

Since dr{x)/dx and s(y) are rational functions, in the above equation the coefficients 
of the logarithmic terms must vanish. Therefore: 



dr{x) 

dx 
dr{x) 



x^y+{{-C3 -l)x^ - cix-l\y = Q (45) 

^^ x^y+({-C2-2)x^ + {2c3 + 2c2 + Q + ci)x-C3-C2-2)y = Q (46) 

Solving ( p5|) and ( |45| ) for dr{x)/dx and integrating in x, one gets: 

r(x) = —X — c5 X + x~ — ci ln(x)) (47) 

c2 c3 
r(x) = i-2c3 -2c2-Q- cl) ln(x) + cSx + 2x-2x"^ (48) 

X X 

Imposing that the logarithmic terms vanish and that the right-hand side of (^) and 
(|4q ) are equal, we get: 

{cl =Q,c2 = -c3 -2,,c3 = c3}. (49) 



Using these in (2£), solving the resulting equation for ds{y)/dy and imposing that 



it can not deppend on x, we have: 

(1 + c3) x^ + ((2 + 2 c5) y - 2 - 2 c5) X + (1 + c3) y"^ + {-2 - 2 c3)y + 1 + c3 = 0, 

(50) 
implying that C3 = —1. So, using (^9|): 



{ci =0, cS = -2, c5 = -1}. (51) 

Substituting (^) into (p7[), we get r(x) = 1/x. Using this in eq.(^) and solving for 
s{y), we have finally: 

-+- 

r{x) = -, s{y) = -, ^ R= ^f ' (52) 

X y y\x + y) 

5 Conclusion 

In this paper, we have presented an extension of the Prelle-Singer method which 
solves a class of LFOODEs missed by the PS-approach. Although our method deals 
with a restricted class of LFOODEs (of the type (||) and among those the ones 
defined by the three cases considered on section |4.2|) and we are still working on the 



demonstration of conjecture 1 (see section 4.1) that would prove the generality of 



the method, we believe the method to be a valid contribuition since we could not 
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find any counterexample where the conjecture fails and, apart from that, it solves 
some LFOODEs that "escape" from many powerfull solvers, using any method of 



solution (the example in section 4.3.3 is such a case). We hope, in the near future, 
to present results concerning the extension of our method to deal with a larger class 
of LFOODEs. For example, considering ro(x, y) to be a general rational function 
or LFOODEs with M{x, y) and/or N{x, y) not restricted to be polynomials, etc. 
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